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Abstract. We have investigated the structure of evolved giant 
stars with masses 3-10 Mq in order to evaluate the bind- 
ing energy of the envelope to the core prior to mass transfer 
in close binary systems. This binding energy is expressed by 
a parameter A which is crucial for determining the outcome of 
binaries evolving through a common envelope (CE) and spiral- 
in phase. We discuss the A-parameter and the efficiency of 
envelope ejection in the CE-phase, and show that A depends 
strongly on the evolutionary stage (i.e. stellar radius) of the 
donor star at the onset of the mass transfer. The existence of 
this relation enables us to introduce a new approach for solving 
the energy equation. For a given observed binary system we 
can derive a unique solution for the original mass and age of 
the donor star, as well as the pre-CE orbital period. 
We find that the value of A is typically between 0.2 and 0.8. 
But in some cases, particularly on the asymptotic giant branch 
of lower-mass stars, it is possible that A > 5. A high value of 
A (rather than assuming a high efficiency parameter, t^ce > 1) 
is sufficient to explain the long final orbital periods observed 
among those binary millisecond pulsars which are believed to 
have evolved through a CE-phase. 

We also present a tabulation of A as a function of stellar radius 
and mass, which is useful for a quick estimation of the orbital 
decay during a common envelope and spiral-in phase. 

Key words: stars: evolution - stars: mass loss - binaries: gen- 
eral - methods: numerical - pulsars: PSR J1454-5846 



1. Introduction 

A very important stage of the evolution of close binaries is the 
formation of a common envelope. Generally, if a star in a binary 
fills its Roche-lobe (either because of the swelling up of the en- 
velope due to exhaustion of nuclear fuel at the centre - or re- 
duction in the separation due to loss of orbital angular momen- 
tum) there will be an attempt to transfer mass to the companion. 
However, if the thermal re-adjustment timescale of the accret- 
ing star is larger than the mass-transfer timescale, the accreted 
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layer piles up above the companion, heats up and expands so 
that the accreting star will also fill its Roche-lobe. The result 
is that the transfered matter will form a common envelope em- 
bedding both stars (Paczynski 1976; Ostriker 1976). When the 
mass ratio is large the CE-phase is accompanied by the cre- 
ation of a drag-force, arising from the motion of the compan- 
ion star through the envelope of the evolved star, which leads 
to dissipation of orbital angular momentum (spiral-in process) 
and deposition of orbital energy in the envelope. Hence, the 
global outcome of a CE-phase is reduction of the binary sepa- 
ration and often ejection of the envelope. For a general review 
on common envelopes, see e.g. Iben & Livio (1993). 

There is clear evidence of orbital shrinkage (as brought 
about by frictional torques in a CE-phase) in some observed 
close binary pulsars and white dwarf binaries (e.g. PSR 
1913H-16, L 870-2 - see Iben & Livio 1993 and references 
therein). In these systems it is clear that the precursor of the 
last-formed degenerate star must have achieved a radius much 
larger than the current orbital separation. Since the present 
short separation of degenerate stars in some binaries can not be 
explained by gravitational wave radiation or a magnetic wind 
as the main source of the loss of Joib, frictional angular mo- 
mentum loss was probably responsible for the very close orbits 
observed. An alternative process recently investigated (King & 
Begelman 1999; Tauris et al. 2000) is a highly Super-Eddington 
mass transfer on a subthermal timescale. However, as shown 
by the latter authors this process can not explain the observed 
properties of all systems. They conclude that CE and spiral-in 
evolution is still needed to account for the formation of double 
neutron star systems and mildly recycled pulsars with a heavy 
white dwarf companion and Porb < 3 days. 

There are many uncertainties involved in calculations of 
the spiral-in phase during the CE evolution. The evolution is 
tidally unstable and the angular momentum transfer, dissipa- 
tion of orbital energy and structural changes of the donor star 
take place on very short timescales (lO'' — 10'' yr). A complete 
study of the problem would also require very detailed multi- 
dimensional hydrodynamical calculations which is beyond the 
scope of this paper. Here we aim to improve the simple for- 
maUsm developed by Webbink (1984) and de Kool (1990) for 
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estimating the orbital evolution of binaries evolving through a 
CE and spiral-in phase. An important quantity in their formal- 
ism is the binding energy between the envelope and the core 
of the donor star. This quantity is expressed by a parameter A 
which depends on the stellar density distribution, and conse- 
quently also on the evolutionary stage of the star. With a de- 
tailed evolutionary calculation, Bisscheroux (1998) found that 
A ^ 0.4 — 0.6 on the red and asymptotic giant branches of a 
5 Mq star In the literature this parameter is usually taken to 
be a constant (e.g. A = 0.5) disregarding the evolutionary stage 
(structure) of the donor star at the onset of the mass transfer 
process. This in turn leads to a misleading (over)estimation of 
the so-called efficiency parameter, rjcE, often discussed by var- 
ious authors. From observations one is able to put constraints 
on the product jycE A (see eq. ^ below). However, by underesti- 
mating A the derived values for t^ce will be overestimated - e.g. 
van den Heuvel (1994) and Tauris (1996) concluded t^ce > 1 
by assuming A = 0.5. By calculating A as a real function of 
stellar radius we are able to improve the estimates of post-CE 
orbital separations and hence put more realistic constraints on 

VCE- 

Detailed calculations on the binding energy of the envelope 
to determine the outcome of binaries evolving through a CE- 
phase have been done by Han et al. (1994, 1995). The method 
and the stellar evolution code used to calculate the binding en- 
ergy in their work is similar to the one applied in our work. 
However, here we discuss the binding energy within the context 
of the A-parameter in the Webbink-formalism. We present cal- 
culated values of A for stars with initial mass in the range of 3 - 
10 Mq at different evolutionary stages. Such values are usefull 
for an estimation of the orbital decay during a common enve- 
lope and spiral-in phase. We also demonstrate a new approach 
to find a unique solution for the stellar and binary parameters 
of systems evolving through a CE and spiral-in phase. 

We will in particular consider binaries in which the com- 
panion of the evolved star is a neutron star. Such binaries are 
the progenitors of binary millisecond pulsars (BMSPs) with a 
massive degenerate companion (e.g. van den Heuvel 1994) and 
they may appear in nature as short lived intermediate-mass X- 
ray binai-ies (IMXBs). 

In Sect. ^ we describe the basic concepts of common en- 
velope evolution concerning the energy equation and the prob- 
lems of defining the binding energy. Sect. ^ introduces briefly 
the numerical stellar code used to evaluate the parameter A. 
The results are presented in Sect. ^ and discussed in Sect. ||. 
We summarize our conclusions in Sect. ^ 

2. Common envelope and spiral-in phase 

2.1. The energy equation and orbital evolution 

A simple estimation on the reduction of the orbit can be found 
by simply equating the difference in orbital energy (before 
and after the CE-phase) to the binding energy of the envelope 
of the (sub)giant donor. Following the formalism of Webbink 
(1984), the original binding energy of the envelope at the onset 



of mass transfer will be of the order: — GAf donor -^^cnv/oi'^L, 
where A/donor is the mass of the donor star at the beginning of 
the CE-phase; M^m is the mass of the hydrogen-rich envelope 
of the donor; Oi is the orbital separation at the onset of the CE 
and tl — Rh/o,i is the dimensionless Roche-lobe radius of the 
donor star, so that airL = i?L ~ ^^donor- A parameter A was 
introduced by de Kool (1990) as a numerical factor (of order 
unity) which depends on the stellar density distribution, such 
that the real binding energy of the envelope can be expressed 
as: 

G A/donor A'/cnv 



E. 



AajrL 



The total change in orbital energy is given by: 

GAfcoreMs , GA^d onor M2 



AEorh = -- 
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where Afcorc — Aidonor — Afenv is the mass of the helium core 
of the evolved donor star; M2 is the mass of the companion star 
and Of is the final orbital separation after the CE-phase. 

Let yycE describe the efficiency of ejecting the envelope, i.e. 
of converting orbital energy into the kinetic energy that pro- 
vides the outward motion of the envelope: -Ecnv = Vce Ai?orb 
or, by equating eqs. (|l]) and (||): 
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2.2. The binding energy of the envelope 

The total binding energy of the envelope to the core is given 
by: 



A/cor. 



iJbind = / ( -^^ + [7 Mm 



(5) 



where the first term is the gravitational binding energy and U 
is the internal thermodynamic energy. The latter involves the 
basic thermal energy for a simple perfect gas (35RT/2/i), the 
energy of radiation (aT*/3p), as well as terms due to ionization 
of atoms and dissociation of molecules and the Fermi energy of 
a degenerate electron gas (Han et al. 1994, 1995). 

The binding energy of the envelope can be expressed by 
eq. (^, under the assumption that the entire internal energy 
is used efficiently in the ejection process. However, it is not 
known how much of the internal energy is involved to unbind 
the envelope. Han et al. (1995) introduced a parameter ath and 
expressed the envelope binding energy as: 



Enrw 



GM{r) 



dm + Qfth 



A/c, 



Afd, 



Udm 



(6) 



A/c, 



The value of ath depends on the details of the ejection process, 
which is very uncertain. A value of ath equal to or 1 cor- 
responds to maximum and minimum envelope binding energy. 
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respectively. By simply equating eqs. and (|]) we calculated 
the parameter A for different evolutionary stages of a given star 
The minimum and maximum derived values of A are denoted 
by Ag (ath = 0) and Ab (ath = !)■ 

By taking the core mass Mcoro as the lower boundary of 
the integral, we assume that the structure of the core does not 
change during the envelope ejection process, and therefore we 
assume that there is no exchange of energy between the core 
and the envelope. Here, we define the core mass as the central 
mass which contains less than 10 % hydrogen - see discussion 
in Sect. IJ. 

For the case of binary stars we use the radii of the stars 
to calculate the separations at the onset of the mass transfer 
assuming a companion star of mass M2. An estimation on the 
Roche-lobe radius is given by Eggleton (1983): 

^ 0.49 gj 

^ 0.6 + g-2/3 1n(l + gi/3) ^ ' 

where q — Mdonor/^^2 is the mass ratio. 



Table 1. The estimated A (see text) as a function of stellar ra- 
dius for different donor stars (graphs shown in Fig. |l]). i? is in 
units of Rq and — I /ME? is the gyration radius of the star 
The corresponding orbital period, Porb (days) is also given at 
the onset of Roche-lobe overflow, assuming a 1.3 Mq neutron 
star companion. 



M = 3.0Mq M = 4.0Mq 

R Ag Ab k^ Porb R Ag Ab k^ Porb 

20.0 * 0.46 0.87 0.13 16.4 20.0 0.17 0.33 0.02 13.6 

40.0 * 0.47 0.94 0.14 46.4 40.0 * 0.34 0.64 0.13 38.6 

70.0 ** 0.46 0.97 0.12 107.6 70.0 " 0.39 0.76 0.12 89.2 

100.0 ** 0.42 0.94 0.12 183.9 100.0 " 0.43 0.91 0.12 152.6 

200.0 ** 0.48 1.81 0.12 519.5 200.0 " 0.59 2.32 0.12 432.1 

300.0 " 0.49 3.44 0.12 953.5 300.0 " 0.64 5.43 0.13 791.0 

400.0 " 0.50 16.50 0.13 1470.0 400.0 " 0.59 10.60 0.14 1220.0 

500.0 " 0.51 0.14 2048.9 500.0 " 0.61 0.15 1703.4 

600.0 " 0.54 0.15 2700.7 600.0 " 0.63 0.16 2242.4 



3. A brief introduction to the numerical computer code 

We used an updated version of the numerical stellar evolution 
code of Eggleton (1971, 1972, 1973). This code uses a self- 
adaptive, non-Lagrangian mesh-spacing which is a function 
of local pressure, temperature, Lagrangian mass and radius. It 
treats both convective and semi-convective mixing as a diffu- 
sion process and finds a simultaneous and implicit solution of 
both the stellar structure equations and the diffusion equations 
for the chemical composition. New improvements are the in- 
clusion of pressure ionization and Coulomb interactions in the 
equation-of-state, and the incorporation of recent opacity ta- 
bles, nuclear reaction rates and neutrino loss rates. The most 
important recent updates of this code are described in Pols et al. 
(1995, 1998) and some are explained in Han et al. (1994). 

We performed such detailed numerical stellar evolution cal- 
culations in our work since they should yield more realistic 
results compared to models based on complete, composite or 
condensed poly tropes. 

4. Results 

The above-mentioned detailed stellar evolution code has been 
used to probe the structure of donor stars at the onset of 
mass transfer. We evolved stars with an initial mass in the 
interval of 3 - 10 Mq assuming a chemical composition of 
{X — 0.70, Z = 0.02) and using a mixing-length parameter 
of a = I /Hp = 2.0. Convective overshooting is taken into 
account in the same way as in Pols et al. (1998) using an over- 
shooting constant So^ = 0.1. We used de Jager's wind mass- 
loss rate (de Jager et al. 1988; Nieuwenhuijzen & de Jager 
1990) to estimate the mass loss prior to the Roche-lobe over- 
flow and CE evolution. The wind mass loss from the donor star 
(and subsequent widening of the orbit) prior to the Roche-lobe 
overflow is typically only ^ 1 — 4 % on the AGB and hence 
this effect is not very significant for the stars considered here. 
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Fig. 1. The parameter A as a function of stellar radius, for stars with initial mass of 3 - 10 Mq. Upper and lower solid lines 
represent the value of A derived from the total (Ab) and gravitational (Ag) binding energy, respectively - see text. For stars with 
mass in the interval of 3 - 6 Mq, the Ab -curv es ar e plotted until the moment before the binding energy becomes positive (and 
the Ab-value would become negative, see Sect. 5.1 for an explanation). The dotted line in each panel presents the gyration radius 
of the star (k^ ~ I/MR^) with a scale given on the right y-axis. The vertical lines separate the regions of case A (left-), case B 
(middle-) and case C (right-side) mass transfer from a binary donor 
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Fig. 2. Total and gravitational binding energy (upper and lower solid line, respectively) around the core of stars with an initial 
mass of 4 Mq (left-side panels) and 10 Mq (right-side panels) at three giant stages. The stellar radius corresponding to each 
stage is indicated at the bottom of the panels. The density profile is represented by a dotted line (with the scale given on the right 
y-axis). Solid and open stars indicate the core mass chosen to determine the binding energy using our choice and the definition 
of Han et al. (1994), respectively. 



With the Eggleton code, the total binding energy is com- 
puted for each mesh-point. Included here are the ionization of 
H+, H, He+"'", He+, He, dissociation of H2, rotational and vi- 
brational modes of H2, as well as ionization of seven heavier 
elements (C, N, O, Ne, Mg, Si, Fe) which are assumed to be 
fully ionized at all temperatures and densities (Pols et al. 1995). 
By integrating from the core to the photosphere, we found the 
total binding energy of the envelope. We used this energy to de- 
rive the parameter A. In Fig. [l| we present the relation between 
A and stellar radius for each of the different donor stars con- 
sidered. Here A was determined either from the total binding 



energy of the envelope (i.e. ath = 1) or from the gravitational 
binding energy alone (i.e. ath = 0)- 

For a donor star in a binary system the stellar radius at the 
onset of the CE-phase is roughly equivalent to its Roche-lobe 
radius, and given a mass of its companion star one finds the 
corresponding orbital period of the system. In Table |l], assum- 
ing that the companion is a 1.3 Mq neutron star, we list the 
relation between the orbital period of the binary and A at the 
onset of the mass transfer (again due to the wind mass loss, this 
period might be slightly larger than the initial ZAMS period). 
We distinguish Ag (as derived from gravitational binding en- 
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ergy alone) from Ab (calculated from the total binding energy). 
The gyration radius of the star, — I/MR^, is also presented 
to provide information about the orbital tidal instability. Mass 
transfer becomes tidally unstable if the rotational angular mo- 
mentum of the donor exceeds one third of the orbital angular 
momentum (Darwin 1908). 



5. Discussion 

5.7. On the binding energy of the envelope 

In Fig. ^ we have plotted the binding energy of the envelope 
as well as the density profile of two stars with an initial mass 
of 4 and 10 Mq. When the star reaches the asymptotic giant 
branch, the core is surrounded by a convective envelope, and 
the density gradient around the core mass is very steep. As the 
star ascends this giant branch, the convective envelope becomes 
deeper and the density contrast higher. 

Assuming a star to be a non-rotating globe of non- 
degenerate ionized hydrogen and helium in hydrostatic equi- 
librium (i.e. nearly an ideal monatomic gas) one would expect 
from the virial theorem: Ab — 2Ag. However at late stellar evo- 
lutionary stages this picture no longer holds. After the second 
dredge-up stage of evolution, the gravitational binding energy 
becomes smaller as the stellar radius expands to a supergiant 
dimension. Meanwhile the ionization and dissociation energy 
becomes more significant when the temperature decreases in 
the envelope (Bisscheroux 1998) - raising the internal energy 
above 0.5 times the gravitational binding energy in absolute 
value. Hence, the total binding energy becomes smaller in ab- 
solute value relative to the gravitational binding energy. This in 
turn causes a larger discrepancy between Ag and Ab, see upper 
right panel in Fig. |l]. 

Especially for low-mass stars (3-6 Mq) we see that the 
contribution of the internal energy to the total binding energy 
becomes dominating. At some point on the asymptotic giant 
branch, the internal energy dominates over the gravitational 
binding energy, and hence the total binding energy of the enve- 
lope reaches a positive value. One might regard this alteration 
from a negative to a positive value of the total binding energy 
as the ejection of the giant's envelope, or as a point of the on- 
set of a pulsationally driven (Mira) superwind (Han et al. 1994 
and references therein). However, when the binding energy be- 
comes positive, the derived parameter Ab will be negative and 
our formalism in Sect. [2.l| breaks down - cf. resulting "nega- 
tive" final separation in eq. (^. On the contrary Ag is relatively 
constant during this late stage of evolution. 



The evolution of stars above 7 Mq is cut short by car- 
bon burning. Therefore they never move as far up the AGB 
as lower mass stars, and the internal (recombination) energy 
never becomes so important as to dominate the total binding 
energy. This is why the difference between Ag and Ab is less 
pronounced as we go to the more massive stars. 




600 



Fig. 3. Ab as a function of radius for a 6 Mq star The full line 
is for X = 0.70, Z = 0.02 and a = 2.0. The dashed line is for 
the same chemical composition with a = 3.0. The dotted line 
is for a star with X = 0.75, Z = 0.001 and a = 2.0. See text 
for a discussion. 



5.2. On the definition of the core mass boundary 

Our results are sensitive to the definition of core mass. In our 
work we have chosen to define the core mass as the central 
mass which contains less than 10 % hydrogen. Using other 
boundary conditions such as log p/dm^ = (e.g. Biss- 
cheroux 1998) or the alternative method by Han et al. (1994) 
will lead to different values for the binding energy of the en- 
velope and hence different estimates of A. In general the above 
mentioned methods will result in a core mass boundary fur- 
ther out in the star compared to our definition (see symbols on 
the dashed line in Fig. ^). Hence their methods lead to smaller 
envelope binding energies and thus larger values of A (factor 
- 2) on the RGB and EAGB. For a 10 Mq star, this difference 
is much higher than a factor of 2. Further up the AGB, where 
the density gradient is steeper, the differences become smaller 
and the result more reliable. Nevertheless, in accordance with 
our findings using any of the two other boundary conditions 
mentioned will also result in large values of A on the AGB (for 
low-mass stars: Ab ^ 1) and this is the important issue. 

5.3. Dependency of the X-paratneter on the chemical 
composition and mixing-length parameter 

In order to investigate the dependency of the A-parameter on 
chemical composition and mixing-length parameter, we have 
plotted in Fig. || the computed Ab (R) for both Pop. I and Pop. II 
chemical abundancies, as well as for two different values of the 
mixing-length parameter for a 6 Mq donor star. 

^From the ZAMS to the base of the RGB (R ^ 50 Mq) 
there are no differences on the calculated values of Ab. It is 
also seen that there is hardly any dependency on the chemical 
composition for the rest of the evolutionary track far up the 
AGB. However, from the horizontal branch (R ~ 150 Mq) 
and halfway up the AGB (R ^ 300 Mq) a higher value of Ab is 
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Fig. 4. Final orbital separations after a CE-phase for systems 
with an initial donor star of mass 9 Mq and 5 Mq (upper and 
lower panel, respectively). The companion is assumed to be a 
1.3 Mq neutron star. The full line is the final separation we 
obtain when the parameter A was calculated as a function of 
stellar radius and the dashed line is for a constant A = 0.5. The 
dotted line is the critical final separation at which the helium 
star (the core of the stripped donor) immediately fills its Roche- 
lobe after the CE and spiral-in phase. If the final separation, at 
is smaller than this value the system is most likely to merge and 
hence the binary will not survive. The radius of the helium star 
was calculated with our numerical stellar evolution code. 



found using a larger value of the mixing-length parameter. This 
picture is reversed near the top of the AGB (R > 400 Mq). In 
general we can conclude that Ab is almost independent of the 
chemical composition of the star, but it depends on the choice 
of the mixing-length parameter. 

5.4. The efficiency parameter of the ejection process 

In order to explain the long orbital periods observed among 
those of the BMSPs which are likely to have evolved through 
a CE-phase, Tauris (1996) and van den Heuvel (1994) found it 
necessary to invoke efficiency parameters larger than unity (e.g. 
1 < '?CE < 4). However, their simple estimates were based on 
a constant value of A = 0.5. Here we suggest that a high value 



of A (rather than a high value of t^ce) is the natural explanation 
for the long orbital periods observed. This in turn indicates that 
if the internal energy is used efficiently (ath = 1) to unbind 
the envelope, there is no need to assume t^ce > 1 for these 
systems. This result was also concluded by Han et al. (1995). 

Let us now consider a system with a 9 Mq donor star and 
a 1 .3 Mq neutron star. Such a binary is the progenitor of a pul- 
sar with a massive degenerate companion. We computed the 
final orbital separation of the binary after a CE-phase by us- 
ing i?onv = f?CE A£^orb Combined with eqs. (||) and (^). Here 
we applied tjce = 1 and = 1 (i-C- including the internal 
thermodynamic energy, A = Ab), and present the result in the 
upper panel of Fig. Q As a comparison, we also show the final 
separation computed by using eq. (Q) with a constant efficiency 
parameter A = 0.5 (commonly used in the literature). 

A binary will survive the CE-phase only if the final orbit is 
large enough so that the core of the donor will not instantly fill 
its Roche-lobe. In the case of a binary with a 9 Mq donor and 
a 1.3 Mq neutron star, this condition is reached if the binary 
separation at the onset of the CE-phase is larger than 435 Rq 
(^oib > 327^). If we use the constant A = 0.5, this separation 
is reduced to 175 Rq (Poib > 84^^). In the latter case the donor 
star still has a radiative envelope with only a very thin convec- 
tive shell at the onset of the mass transfer. Therefore the CE is 
formed only as a result of orbital shrinkage due to mass transfer 
with a high mass ratio between the donor and the neutron star. 
In the former case, the donor star had time to develop a deep 
convective envelope and it will therefore expand in response to 
mass loss due to its isentropic entropy profile. This expansion, 
in combination with the orbital shrinkage mentioned above, 
will immediately lead to a runaway mass transfer, and hence 
a CE will be formed easily. 

At late stages of stellar evolution, the value of Ab has an im- 
portant effect on the final orbital period. This effect is more 
significant in binaries with less massive donors as can be seen 
in the lower panel of Fig. Q for the case of a 5 Mq donor. 

5.5. A new approach for tracking the evolutionary history of a 
system evolving through a CE 

As we have demonstrated in this paper, A is a function of stellar 
radius for any given mass. Combining this fact with eq. (^, 
which relates the observed orbital period to the original radius 
of the Roche-lobe filling donor star as a function of mass and A 
(or equivalently: Of /as), we are able to place severe constraints 
on the mass and radius (and hence age) of the original donor. 

An example of this approach is shown in Fig. ^ where 
we have used the parameters for the recently discovered bi- 
nary pulsar PSR J 1454-5846 (Camilo et al. 2000). This system 
has an orbital period of Porb = 12^42 and a mass function, 
/ — 0.130 Mq. The exact mass of the white dwarf in this 
system is unknown; but as an example we shall assume an or- 
bital inclination angle, i = 60° yielding A/wd = 1.03 MQ(for 
A^NS = 1-30 Mq). Hence we must require that the original 
giant donor star had a core mass of this value by the time its 
envelope was ejected in the CE-phase. In Fig. |^ the intersection 
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Fig. 5. The parameter A = Ab as a function of stellar radius, for donor stars with initial mass of 5 - 7 Mq. Also shown are 
curves for eq. (Q) calculated with the observed parameters given for PSR J 1454-5846 and assuming rycE = 10. The intersection 
of these curves yields a unique solution for the properties of the original donor star (the progenitor of the white dwarf) and the 
pre-CE orbital period - see text for a discussion. 



points for a 5 Mq and 6 Mq donor yield core mass values of 
0.93 Mq and 1.10 Mq, respectively. To obtain a core mass of 
1.03 Mq we must therefore search for a solution in-between. 
We find Mdonor = 5.6 Mq and the radius at the onset of the 
RLO is 340 Rq. The corresponding age of the donor is 93 Myr 
and the pre-CE orbital period is 760*^. 

The above solution was found for an inclination angle, 
i = 60°. In Fig. ^ we have plotted the solutions for Mdonor and 
i?L as a function of inclination angle for two different values 
of the efficiency parameter, 77CE. It is interesting to notice that 
Mdonor is more or less independent of rycE whereas i?L gets 
larger for smaller values of yycE- This result is expected theo- 
retically since a low value of 77CE requires either more orbital 
energy released and/or less binding energy of the envelope, and 
thus a larger i?L, in order to eject the envelope successfully. An 
upper limit to the mass of the white dwarf is about 1.3 Mq 
which yields i > 47°. For 7/ce — 0.5 no solution is found for 
i < 59°. Hence it would be interesting if a future detection of a 
spectroscopic line from the white dwarf would yield an orbital 
inclination angle less than this value, since small values of tiqe 
would be excluded that way. 



It should also be noted from Fig. ^ that the binary pulsar 
must have had an AGB radius at the onset of the RLO. Many 
of the observed binary pulsars which have survived a common 
envelope evolution also have long orbital periods. Hence, the 
density gradient was fairly steep at the onset of the Roche- 
lobe overflow and therefore the uncertainties regarding the core 
mass definition as already discussed in Sect. 5.2 are less impor- 
tant in this context. 



6. Conclusions 

We have adapted a numerical computer code to study the stellar 
structure of stars with an initial mass of 3 - 10 Mq in order to 
evaluate the binding energy of the envelope to the core which 
determines the parameter A in the energy equation of the com- 
mon envelope evolution. We have presented evidence that the 
parameter A depends strongly on the evolutionary stage. Hence 
the value of A, at the onset of the mass transfer and CE-phase, 
is unique for binaries with different initial orbital periods and 
donor masses. Taking advantage of this, we have demonstrated 
a new approach for finding a unique solution to the energy 
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Fig. 6. Solutions for the mass (top panel) and radius (bottom 
panel) of the giant progenitor of the massive white dwarf ob- 
served in PSR J 1454-5846, as a function of orbital inclination 
angle, i and CE efficiency parameter, ijce- On the top is shown 
MwD as a function of the inclination angle - see text for further 
details. 
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equation and the original mass and radius of the donor star. 
An application of this approach on PSR J1454-5846 yields a 
constraint for the donor star and the pre-CE orbital period as a 
function of the unknown orbital inclination angle. 

We have demonstrated that in order to obtain a long final 
orbital period after a CE-phase (as observed in some BMSPs), 
we do not require an efficiency parameter larger than unity, if 
the internal energy can be used efficiently to eject the envelope. 

In order to obtain the final orbital separation of a binary 
system after a CE and spiral-in phase, we advise one to calcu- 
late the binding energy from the stellar structure by means of 
eq. (^ with < ath < 1. In case the detailed structure of the 
donor is not available (e.g. for a quick back-of-the-envelope- 
calculation), eq. (Q) can be used with rjcE = 1 in combination 
with Ag < A < Ab. The value of A can be found in our Table |l] 
for a given donor star mass and radius. 



